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ABSTRACT 


X 


The  blast  wave  solutions  for  the  axlsymmetric  flow  fields 
associated  with  a  family  of  power -law  shocks,  “ - ~n 


modified  to  properly  account  for  the  entropy 
layer  near  the  surface  of  the  body.  The  method  of  inner 
and  outer  expansions  was  used.  It  was  shown  that  thia 
method  is  of  limited  usefulness  for  n  only  slightly  greater 

than  1/2  due  to  the  slow  convergence  of  the  inner  expan - 

*  * 

sions,  but  for  (y'i  2)/(3 X  +  2)5n  «1  the  expansions  con¬ 
verge  rapidly  enough  for  the  method  to  be  useful.  The 
case  n  =  1/2  is  singular  and  causes  no  difficulty. 

In  the  range  Of  +  l)/(2y  +  l)Sn<l  the  body  shape  given 
by  blast  wave  theory  is  correct,  but  the  blast  wave  flow 

field  profiles  must  be  modified  near  the  body.  For 

-A* 

n  «e  (y  +  l)/(2y  +  l)both  the  body  contour  and  the  flow  field 
profiles  must  be  modified.  \ 
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1.  INTRODUCTION 


The  inviscid  hypersonic  flow  of  a  perfect  gas  about  slender  blunt 
nose  bodies  has  been  the  subject  of  numerous  investigations.  Blast  wave 
theory  has  proven  to  be  a  useful  method  for  analysing  this  problem.  This 
approach  uses  the  fact  that,  subject  to  the  hypersonic,  small  disturbance 
conditions,  the  flow  field  about  a  blunt  nose  slender  body  is  similar  to  the 
unsteady  flow  field  of  an  intense  blast  wave  propagating  into  a  medium  at 
rest.  Lees  and  Kubota  (Reference  1)  and  Kubota  (Reference  2)  have  used 
the  blast  wave  analogy  to  analyse  the  inverse  problem  of  a  family  of  power- 
law  shocks,  y#~xn.  Solutions  were  obtained  for  1/2  £n  £1.  A  para¬ 
boloidal  shock,  n  =  1/2,  was  found  to  give  a  sero  thickness  afterbody 
with  all  of  the  drag  concentrated  at  the  blunt  nose.  These  solutions,  how¬ 
ever,  do  not  properly  account  for  the  layer  of  high  entropy  air  surrounding 
the  body.  This  entropy  layer  consists  of  streamlines  that  pass  through  the 
blunt  portion  of  the  shock.  Singularities  such  as  infinite  temperature  and 
sero  density  at  the  body  surface  were  obtained. 

Sychev  (Reference  3)  introduced  the  idea  of  obtaining  a  corrected 
body  shape  which  accounts  for  the  presence  of  the  entropy  layer  and  pro¬ 
duces  the  outer  flow  field  and  shock  wave  given  by  blast  wave  theory.  He 
obtained  a  numerical  solution  for  the  body  that  produces  a  paraboloidal 
shock  and  showed  that  a  critical  value  of  n  exists  above  which  the  body 
shape  is  not  affected  by  the  entropy  layer. 

Yakura  (Reference  4)  has  applied  the  method  of  inner  and  outer  ex¬ 
pansions  as  developed  by  Lagerstrom,  Kaplun,  and  Cole  (References  5 
and  6)  for  treating  singular  perturbation  problems  to  the  inverse  problem 
with  specified  shock  shapes.  He  obtained  solutions  for  blunted  wedges, 
blunted  cones,  and  a  paraboloidal  shock.  The  paraboloidal  shock  gives  a 
body  which  grows  as  x*^  in  the  downstream  direction  instead  of  the 
sero  thickness  afterbody  predicted  by  blast  wave  theory. 

The  present  analysis  uses  the  methods  developed  by  Yakura  to  add 
an  entropy  layer  to  the  blast  wave  solutions  obtained  by  Kubota  for  a  family 
of  power -law  shocks.  Only  the  axisymmetric  case  is  considered,  but  the 
methods  are  applicable  to  plane  flow. 
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2.  METHOD  OF  ANALYSIS 


2.  1  Equation!  and  Boundary  Condition 

In  considering  the  inverse  problem  with  the  shock  wave  specified 
it  is  convenient  to  use  the  Von  Mises  transformation.  The  indeper.dent 
variables  are  x.  a  downstream  coordinate,  and  vp,  a  variable  that  is  pro¬ 
portional  to  the  streamfunction.  These  variables  make  it  easy  to  identify 
the  body  as  the  ip  =  0  streamline,  and  the  conservation  of  entropy  along 
streamlines  helps  to  define  the  entropy  layer. 

Following  the  analysis  of  Yakura,  Reference  4,  the  streamfunction 
for  axisymmetric  flow  is  defined  by: 


0*  0* 

87  =  PUV‘  8x  =  -pVy 


(2.1) 


Introducing  the  Von  Mises  variables  x  and  «p  =  k¥,  where  k  is  an 
arbitrary  constant  that  can  be  chosen  for  convenience,  the  hypersonic 
equations  of  motion  for  a  perfect  gas  may  be  written  as: 


>13 

II 

£i<g 

(2.2) 

•< 

ii 

(2.3) 

Z  ®E  =  o 

8x  k  8<p 

(2.4) 

+  v2  +  E  =  i 

y  -  1  p 

(2.5) 

=  2G  sin2  0 

(2.6) 
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The  coordinate  ayatem  is  shown  in  Figure  1.  In  Equations  (2.  1) 
through  (2.6)  the  variablea  are  nondimenaionalised  as  follows: 


P 


=  JLl 

^oo 


U,  V  = 


CD 


P  U 
r00  CD 


x’.y1 

x>  y  s  ~ 


(2.7) 


Figure  1.  Coordinate  System 


The  parameter  G  is  a  function  of  y, 


(2.8) 


8  is  the  angle  at  which  a  streamline  crosses  the  shock,  and  L  is  a  charac¬ 
teristic  dimension  of  the  flow  field.  The  energy  equation,  Equation  (2.  5), 
has  been  written  in  integrated  form.  The  axial  momentum  equation, 


u 


I  9p  yv  9p 
p  8x  k  8v^ 


0 


(2.9) 


has  been  replaced  by  Equation  (2.6)  which  states  the  condition  of  constant 
entropy  along  streamlines.  The  usual  hypersonic  approximation 

(M^e)2  =-=»  1  0  -=-=l  (2.10) 


is  made  throughout  this  report  such  that  the  solution  is  independent  of 
free -stream  Mach  number. 
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The  boundary  conditions  at  the  shock  may  be  written  as: 


p  = — r— r  sin  0 
*a  v  +  1 


V  +  1 

Pe  =y^T 


UB  =  (y”tt)  ®  +  c°8^  ® 


(2.  ID 


v  = 


TT 


sin  0  cos  0 


yB  =  yB(x) 

where  the  shock  shape,  yB(x),  is  a  specified  function  for  the  inverse  problem. 
We  consider  the  family  of  power -law  shocks 


y  1  =  Ax' 

'  s 


(2.  12) 


where  A  is  an  arbitrary  constant.  The  remainder  of  this  analysis  there¬ 
fore  reduces  to  that  of  Yakura  (Reference  4)  for  the  spacial  case  n  =  1/2. 

A  characteristic  dimension  is  defined  by 

(A  W<l*n> 

L  .  ($.)  «•  «» 

such  that  the  equation  fo*  the  family  of  shocks  is 

yB=V2xn  (2.14) 

For  n  =  1/2,  L  is  the  nose  radius  of  the  shock,  but  the  nose  radius  is  sero 
for  1/2  <n  «1.  However,  the  physical  significance  of  L  for  1/2  £n*  1 
can  be  seen  by  noting  that,  for  y#  =  1 


dys 

~0(1)  and  x~0(l)  (2.15) 
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Thus  L  is  a  measure  of  the  radius  of  the  strong  blunt  portion  of  the  shock 
in. the  nose  region.  The  ~\fz  term  is  included  only  for  convenience  in 
matching  Yakura's  result  for  n  =  1/2  (Reference  4). 

Continuity  considerations  require  that: 


(2.  16) 


Using  this  relation  and  taking 


y  n/(l -n) 

k  =  (4n  )  (2. 17) 

for  convenience,  the  condition  of  constant  entropy  along  streamlines, 
Equation  (2.6),  can  be  expressed  as: 


jg_  _  2G 

pv  '  i  + 


(2.  18) 


Equation  (2.  18)  illustrates  the  convenience  of  using  the  Von  Mises  trans¬ 
formation  for  the  inverse  problem  since  the  entropy  function  is  known  in 
terms  of  the  dependent  variable  ij>.  The  boundary  conditions,  Equations 
(2.  11),  can  be  expressed  as: 


^  fa  _ a _ 

P«  ~ Y  +  *  i  +  2^(1 -n)/n 


v  +  1 

p.  =J-n- 


u  =* 


-  1 


Z\ \i 


( l-n)/n 


8  =  Y  +  1  l  +  +  1  +  24,<1'n>/n 


(2.  19) 


.  *  K1-,,n) 

VB  7+"1  1  +  2^1-nJ,n 


V,  =  (2k+a) 


1/2 
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The  problem  ia  now  specified  by  Equations  (2.2),  (2.3),  (2.4),  (2.5), 

(2.  18)  and  the  above  boundary  conditions  for  p,  p,  u,  v,  y  as  functions 
of  x  and  4>. 

2.  2  Method  of  Inner  and  Outer  Expansions 

The  method  of  inner  and  outer  expansions  was  developed  by 
Lagerstrom  et.al.,  (References  5  and  6)  for  treating  singular  perturba* 
tion  problems  with  particular  emphasis  on  viscous  flow  at  high  and  low 
Reynolds  number.  Van  Dyke  (Reference  7)  used  the  method  to  identify 
certain  second  order  boundary  layer  effects  such  as  curvature,  vorticity 
interaction,  etc.  Yakura  (Reference  4)  has  used  the  method  to  analyse 
the  inviscid  hypersonic  flow  about  slender  blunt  nose  bodies  accounting 
for  the  entropy  layer  near  the  body.  Essentially,  the  method  divides  the 
flow  field  into  two  regions  in  which  different  asymptotic  expansions  are 
valid  and  uses  a  limiting  process  to  match  the  expansions  in  an  overlap 
region. 

The  method  as  applicable  to  the  inviscid,  blunt  nose,  slender  body 
problem  is  described  below  in  a  way  that  permits  an  understanding  of  the 
present  analysis.  A  more  general  and  detailed  description  may  be  found 
in  the  above  references. 

The  outer  region  of  thq  flow  field  is  defined  as  consisting  of  those 

I 

streamlines  which  cross  the  weak  portion  of  the  shock  where  x»*l.  In 
this  region  i|<  ~0(v)jg)  '~0(x2n).  An  independent  variable  of  order  one  is 
defined  for  the  outer  region  by 

w=~fe~°(i)  U.20) 

x 

and  an  outer  expansion  for  the  dependent  variables  is  assumed  of  the 
form 


Oex  g(x,  «|<)  =  5]  E.(x)F.(w)  (2.21) 

i=l  1  1 

The  functions  E^(x)  can  be  chosen  to  match  the  known  shock  conditions 
expanded  for  large  x,  and  the  boundary  conditions  on  F^(w)  are  taken  at 
the  shock.  This  expansion  scheme  leads  to  an  outer  solution  that  is 
asymptotically  valid  for  large  distances  downstream  from  the  nose. 


Streamlines  near  the  body  cross  the  strong  portion  of  the  shock 
near  the  blunt  nose  and  define  an  inner  region  where  —0(1).  The  ex¬ 
pansion  assumed  above  where  E^(x)  is  taken  from  the  shock  conditions 
for  large  x  is  not  valid.  An  expansion  of  the  following  form  is  assumed 
for  the  inner  region: 

lex  g(x,<)<)  =  Y,  e.(x)f,(«|»)  (2.22) 

1=1  1  1 

Since  the  body  shape  is  unknown  a  priori,  the  form  of  e^(x)  and  the 
boundary  conditions  on  f^(tj»)  must  be  taken  from  the  asymptotic  behavior 
of  the  outer  solution  as  the  inner  region  is  approached.  Two  limit  proc¬ 
esses  are  defined: 

The  Inner  Limit  of  the  Outer  Expansion 

Ilim  Oex  g(x,  t|<)  =  lim  [Oex  g  (x,  4>fJ  (2. 23) 

x— *oo,  4<  fixed 

The  Outer  Limit  of  the  Inner  Expansion 

Olim  lex  g(x,  vJj)  =  lim  [lex  g  (x,  ijj)]  (2. 24) 

x  — ►  oo,  w  fixed 

If  the  outer  solution  is  known,  the  e^(x)  terms  may  be  determined  by  taking 
the  inner  limit  of  the  outer  expansion.  The  assumed  inner  expansion, 
Equation  (2.  22),  may  then  be  used,  and  the  boundary  conditions  on  fj(i|») 
are  taken  such  that: 

Olim  lex  g(x,  4»)  =  Ilim  Oex  g(x,  i|i)  (2.  25) 

The  outer  and  inner  solutions  can  be  combined  in  such  a  way  as  to 
provide  a  composite  expansion  that  is  uniformly  valid  across  the  entire 
shock  layer.  Essentially  it  is  desired  to  eliminate  the  incorrect  portions 
of  both  the  outer  and  inner  expansions.  This  can  be  done  in  the  following 
manner: 


g(x,4»)  «  Oex  g(x,t|/)  +  lex  g(x,  i|«)  -  Olim  lex  g(x,i|i)  (2.26) 

Note  that,  due  to  the  boundary  condition  imposed  upon  the  inner  solution, 
the  outer  and  inner  limits  of  the  composite  expansion  reduce  to  the  outer 
and  inner  expansions  respectively. 


-7- 


3.  SOLUTIONS  FOR  POWER -LAW  SHOCKS 


C 


3.  1  Outer  Solution 

2n 

Introducing  the  variable  w  =  iJj / x  into  Equations  (2.2)  through 
(2.  5),  (2.  16),  and  (2.  19)  and  expanding  for  large  x,  the  equations  and 
boundary  conditions  for  the  outer  region  can  be  expressed  as  follows: 


ey 

8x 


2n 


u  8y  _  v 


8co 


u 


(3.1) 


8 v 
8x 


“k 


8y  _  k 
8u  pu 


.  2n£  +  -JU  * 


0 


(3.2) 

(3.3) 


u2  +  v2  +  2^  i  E-  =  1 

y  -  1  p 


(3.4) 


f  "zJ'Wn* 

2^wx  ) 


with  boundary  conditions  at  the  shock,  wg  =  1/k; 


w  +  »(“/") 


(T-nlTn 


tn5T7K  +  " 


(3.5) 


y  +  1 

p. =  y^~t 


u  =  1  -■ 


s  y  +  1 


(3.6) 


VL 


(Y  +  l)(wgx2nj 


(I"ny72n 


1  - 


( 1 -n)/n  + 


yg  * (2kw#)1^2  xn 
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The  shock  conditions  suggest  the  following  form  for  the  outer 
expansions: 

p  =  PjMx'2*1'"*  +  P2(«)x"4(1'n)  +  ••• 
p  =  D^w)  +  D2(«)x‘2^'n^ 

(3.7) 

u  =  1  +  Ul(«s)x“2<l"n>  +  U2(w)x'4(l*n)  +  ••• 

v  =  V1(w)x‘1(1‘n)+  V2(o»)x‘3(1‘n)  +  •  •  • 
y  =  Y  j(w)  x11  +  Y2(w)x '^2  3n^  +  •  •  • 


Substituting  these  expansions  into  the  outer  equations,  Equations  (3. 1) 
through  (3.5),  the  following  set  of  first  order  outer  equations  and  boundary 

conditions  are  obtained: 

dYl  V1 

Yl  *  ^  ST  s  — 

(3.8) 

»  dYl  k 

1 1  dw  ~ 

(3.9) 

,  dV,  dP, 

vi  +~ — jz~  w  sr  Yi^r  *° 

(1  -  n)  k 

(3.10) 

P1  G 

w(i'n,/n 

(3.11) 

2  2v  P 1 

2Ui +  vi  8  0 

(3.12) 
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with  boundary  conditions  at  the  shock,  u»g  =  1/k: 


l  7TT 


.  v  £  1 

h 


u 


s  .Jsl 

Is  y  + 


<2i3/2n 

V.  *  Mi— 

1  Y  +  1 


Y,  « 


Y 

V5- 


(3. 13> 


Note  that  the  first  four  of  the  above  equations  do  not  contain  the  axial 
velocity,  U^,  and  can  therefore  be  solved  independently  of  the  energy 
equation,  Equation  (3. 12).  This  corresponds  to  the  blast  wave  approxi¬ 
mation  of  u  <*»uqo.  The  axial  velocity  can  be  obtained  from  Equation  (3. 12) 
after  the  first  four  equations  are  solved.  The  solution  of  Equations  (3.8) 
through  (3. 11)  can  be  taken  directly  from  Kubota's  blast  wave  analysis 
(Reference  2)  in  terms  of  the  similarity  variable  a  *  y/y#t 


=  2n2F(s)  (3. 14) 

Dj-tfe)  (3.15) 

Vl=VW<s)  (3.16) 

»  Vz  s  (3.17) 


The  functions  F(s),  p(s),  y(e)  are  given  in  numerical  form  in  Reference  2 
for  y  *  1.4  and  n  ■  1/2,  4/7,  2/3,  3/4,  and  1.  An  approximate  analytic 
solution  for  these  functions  is  given  by  Mirels  (Reference  8). 
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Finally,  the  first  order  outer  solution  or  the  blast  wave  solution 
is  given  by: 


p  =  2n2F(s)  x”2(1'n) 

(3.  IS) 

p  =  p(») 

(3.19) 

v  =  V ifn*U )  x"(l"n) 

(3.20) 

(3.21) 

y  =  yi  *xn 

(3.22) 

The  above  results  are  not  valid  near  the  blunt  nose  or  vary  far  downstream 

where  the  hypersonic,  small  disturbance  conditions  are  not  satisfied. 

They  are  also  invalid  near  the  body  and  give  zero  density  and  infinite 

temperature  on  the  surface  of  the  body.  Also,  for  n  *  1/2,  the  body  is 

of  zero  thickness.  These  last  two  difficulties  can  be  eliminated  by  using 

* 

a  different  expansion  near  the  body -in  the  inner  region  or  entropy  layer. 

3.2  Inner  Solution 

The  outer  expansion  assumes  i|i»lp  and  it  can  be  seen  from 
Equation  (2.  18)  that  this  leads  to  a  singularity  at  the  body  surface.  Dif¬ 
ferent  expansions  must  be  used  for  the  inner  region  where  tp  ~0(1).  The 
expansion  of  Equation  (2.  22)  is  valid  if  the  e^x)  terms  are  determined 
by  taking  the  inner  limit  of  the  outer  expansion  according  to  Equation  (2.23). 
Using  this  expansion,  a  solution  of  Equations  (2.  2)  through  (2.  5)  and  (2. 18) 

can  be  obtained  in  terms  of  the  variables  x  and  ip> 

I 

To  obtain  the  Him  Oex  g(x,  ip)  it  is  necessary  to  know  the  outer  solu¬ 
tion  in  analytic  form  near, the  body.  This  can  be  obtained  by  expanding 
the  outer  solution  for  the  pressure  in  a  Taylor  series  about  the  body, 
u>  =0,  and  using  the  resulting  pressure  distribution  to  solve  the  remaining 
outer  equations.  Thus: 

p  »  |pj(0)  +  ("3^^  w  +  •  •  •  J  x-2<l-n>  +  P^w)  *  ^  +  . . .  (3.23) 
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I 


Since  4<  ~ 
becomes: 

P 


0(1)  in  the  inner  region,  o>  —0(1 


x'2  + 


/x2n)  and  Equation  (3.23) 
••  +  P2(u)x'4(l'n)  +  ••• 


(3.24) 


where  (dP^/du)Q  is  finite  and  of  order  one.  From  Equation  (3.24)  it  is 
seen  that  to  the  accuracy  involved  in  the  first  order  solution  the  pressure 
may  be  taken  as  constant  across  the  inner  region  and  equal  to  the  value 
given  by  the  outer  solution  at  the  body  surface.  Substituting 


P^u)  =  Pj(0) 


(3.25) 


into  the  first  order  outer  equations  gives  the  following  results  valid  near 
the  body: 


-n)/ny 


YjM  = 


(3.26) 


1/2 


3.27) 


(3.28) 


(3.29) 


Using  the  above  relations,  the  inner  limit  process,  Equation  (2.23), 
can  be  applied  to  the  outer  solution.  The  result  in  terms  of  the  inner 
variables  can  be  written  as: 


Him  Oex  p(x,  i|*)  =  P^(0)  x  2^ 


(3.30) 


( 
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Him  Oex  p(x,  t|i) 


Ilim  Oex  y(x,  ip)  = 


/  2kv  o  'I1'*  i7('tl-K)JiaTs!  ,  T 

<— I  -  I7I|F^5lJ  *  > 


(3.31) 

.1/2 


(0)2x2n 


(3.32) 


Similar  expressions  can  be  obtained  for  v  and  u. 


An  important  feature  of  the  application  of  the  method  of  inner  and 
outer  expansions  to  this  problem  can  be  seen  from  Equation  (3.  32). 

For  n  =  1/2,  Yj(0)  =  0  (i.  e.,  the  body  given  by  the  outer  solution  has  aero 
thickness),  and  the  Ilim  Oex  y(x,ip)  reduces  to: 


Ilim  Oex  y(x,  ip)  = 


(3.33) 


Equation  (3.  33)  gives  the  following  form  for  the  inner  expansion  which 
is  the  same  as  that  used  by  Yakura  (Reference  2): 


lex  y(x,  >p)  =  Yj(«|i)  x1/2y 


(3.34) 


For  n  =»l/2,  however,  Y^(0)  ^  0,  and  Equation  (3.  32)  must  be  expanded 
resulting  in: 

Ilim  Oex  y(x,  «p)  =  Y^O)  x” 


+ 


f  -3n+ 
+  0  [x 


4(  1 -n)~ 

Y  J  + 


+  Y2(w)x 


-(2-3n) 


+  •  •  • 


(3.35) 


where  the  term  Y2(u>)x  ^2  ^  represents  the  contribution  of  the  second 
order  outer  solution.  Equation  (3.  35)  gives  the  following  form  for  the 
inner  expansion: 
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lex  y(x,  ij/)  =  yj(4»)  x  +  y2(t|<)x  Y  +  y3(tji)x  Y 

+  •••  +  o[x'(2"3n,]+  •••  (3.36) 

To  obtain  an  inner  solution  that  is  consistent  with  the  first  order  outer 
solution,  enough  terms  must  be  retained  in  Equation  (3.36)  that  the  order 
of  the  dropped  terms  is  less  than  or  equal  to  x  ^  3n^.  It  can  be  seen 
that,  as  n-sl/2,  an  increasing  number  of  terms  must  be  retained  in  the 
first  order  inner  expansion.  Also,  for  n  =  1/2,  Yj(0)  =  0  and  Equation 
(3.  35)  is  not  valid.  Thus  Equation  (3.  36)  can  be  used  for  all  n  •* 1/2  if 
enough  terms  are  retained  in  the  expansion. 

The  following  forms  of  the  inner  expansions  are  consistent  with  the 
first  order  outer  solution  in  the  range  of  values  of  n  indicated.  The 
n  »  1/2  case  was  solved  previously  by  Yakura  (Reference  2)  and  is  in* 
eluded  here  only  for  the  sake  of  completeness. 

j  Sn  <1  p(x,  t)»)  »  pj(4i)  x"2^1  *n* 

-2(1 -n) 

p(x,  4<)  =  pj(ij<)  X  Y  (3.37) 

-2(l-n)^i 

u(x,  4<)  =  1  +  Uj(l|<)  X  Y 


V  +  1 

T7TT 


—  n 


y(x,+)  =  yj(+)  xn 
v(x,  >);)  =  Vj(i)j)  x’(1'n) 


(3.36) 


y  *  ^  £  n  *  _ i  ,  v  __n  .  __  ,,i  „-n+ 


37TT 


.  _  2(  1  -n) 

2v  +  1  y(x,t|;)  =  y  1(4»)  *n  +  y^)  x*n  Y 


(3.39) 


v(x,  *J»)  =  Vj(»|»)  x  ^  +  v2(4<)  X  '  ^  ^ ' 
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1 

n  =  I 


y(x,4/)  =  y  x 


(3.40) 


1/(2y) 


v(x,4l)  =  Vj(4l)  X 

For  1/2  c  n  «  (y  +2 )/(3y  +  2)  additional  terms  would  be  required  in  the 
expansions  for  y  and  v.  Only  the  above  intervals  are  considered  in  the 
remainder  of  this  report. 

Substituting  the  above  expansions  into  the  complete  hypersonic 
equations,  Equations  (2.  2)  through  (2.  5)  and  (2. 18),  results  in  the  follow 
ing  system  of  first  order  inner  equations: 


n  1 


1 


=  0 


dp 

"3ijr 

dyj 

W  =  0 

,  -ny,  =  0 


2G 


1  +  24 » 


( 1 -n)/n 


(3.41) 

(3.42) 

(3.43) 

(3.44) 


u 


1 


(3.45) 


n 


V  +  \ 
2y  4  1 


Two  additional  equations  are  required: 


dy2  k 
~W  '  P|Y| 


(3.46) 


(3.47) 
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n  =  1  /2 


Equations  (3.41),  (3.44),  and  (3.45),  plus  the  following  two 


equations  are  required: 


dVl  _  1 

dll'  “  PjYj 

y i =  2^vi 

The  boundary  conditions  are  such  that: 

Ilim  Oex  g(x,  t|<)  =  Olim  lex  g(x,  ij<)  (3.50) 

The  above  system  of  equations  are  easily  solved  and  the  first  order 
inner  solution  can  be  written  as: 


(3.48) 

(3.49) 


<‘ 


p(x,  i)j)  =  P  ^(0)  x 
|1/V 


-2(1 -n) 


(3.51) 


p(x,w.[^]  ‘-I* 


(3.52) 


u(x,i|<)  =  1  -  \  t  Pt(0) 


r  2g  ~i1/v _ 1  x 

L^J 


-2(1 -n)*^- 


V  +  1  —  1 

irn  -  n  *  1 


y(x, ij/)  =  Y  j(0)  x 
v(x,  «|i)  =  nY j(0)  x'(1'n) 


(3.53) 

(3.54) 

(3.55) 


k 


+  2  _  Y  +1 

T7  -  n  *  — 


2y  +  1  y(x,  \\i)  =  Y ^(0)  x  +  y2(4»)  X 


-n+^A 

Y 


(3.56) 


v(x,i|i)  =  nY j(0)  x‘(1'n)  + 


(3,57) 
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where 


( 


■  7^(pf&]  «^'  +  Yi(0|(7+ ,.!)[?>]  *' 


(3.58) 


I(^;n,v)  =  f 
JiL 


co 


r  . 

l 

c 

'e 

i 

c^- 

«M 

.."rTj 

(l  *  2+"-»)'")l,V_ 

di|i  (3.59) 


n  =  2 


(  \l/2r  -|2V 

y(x.4<)=(yVrj  [jv  +  TiP\{6)J  (2'),+  1) 


v-1  1 

+  l)Vx27 


(3.60) 


v(x,  4»)  = 


2[y(v  +  1>]1/Z 

3.  3  Composite  Expansions 


2  12Y 

(y  +  Op 


^Tol  (2*  + 


&-H) 


1)‘T  X 


(3.61) 


Equation  (2.26)  can  now  be  applied  to  the  inner  and  outer  solution* 
resulting  in  a  solution  that  is  uniformly  valid  from  the  shock  to  the  body. 
Thus  the  final  results  are: 


i  = »  «  1 


.2»2FWx-2"-> 


»  ■  +[^er]  ♦  2*  ”  )  -  ° 

u  = 1  -  §$  +  x 
p,w) 


1  l 

1-nJ  ,  "1  '2(l-n) 


(3.62) 


(3.63) 


X  *2.  .  -2(1 -n) 


r  2o  “ 

i/v 

1 

1 

^(i-nj/n 

17 y 

'  '  Y 

(3.64) 


{ 
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i: 


Y  +  1  . 

WT1  -  n  "  1 


y  =  yi  zxn 
a[z  n^(  ss)  x  ^ 


(3.65) 

(3.66) 


Ml-n) 


^rr2  -  n  *ci±TT  y  =  V2  +  [p^rr]  Y 

- 1 -  1  L  1  J  (3.67) 


_/r  ...  -(1  -n)  .  ["2(1  -  n)  1  k  fzc  IV^K-J 

V  =  V2  n0(z'  x  +  — r -  •  n  y  TO)  P  (6)  M'v.n.Y'  x 

L  J  1  L  1  J  (3.68) 


n  =  2  y 


=  V2  zx1/2  +  (VVt)  (7"+  TfP]7o) 


12y 


V1  ^ 

(2*  +  1)2y  -  (2+)^ 


1 

.2? 


0(z)  -1/2  . 

v  =  ~  x  +  — 

yi  2[V(y  +  d] 


1_ 

hyp77  [iv  +  iff,!")]  h*1’ 


■Ui  v-l 

(1|V-U(IV 


X 

(3.69) 

,-H) 

(3.70) 


In  the  above  equations  the  relation  between  the  similarity  variable  z  and 
the  Von  Mises  variable  is  given  by: 


s|3U) 


[z  -  <t>(z)  j: 


2n 


(3.71) 


L 
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4.  NUMERICAL  RESULTS 


The  body  shape h  associated  with  the  family  of  power-law  shocks 
can  be  obtained  by  setting  ^  ^  in  the  inner  solutions,  Equations  (3.51) 

through  (3.61).  The  results  are: 


V  ±  1  «= 

2y  +  1  ~ 


n 


1 


=  Y^O)  xn 


(4.1) 


V  4  l  <  n  «  -X-Ll 
3y  ±  2  ~  2y  +  1 


yb  =  Yj(0)xn  + 


2G 

L^J 


l/y 


-n+ 


2d -n) 


l(0;n,y))x 


(4.2) 


1 


(4.3) 


The  effect  of  the  entropy  layer  on  the  body  shape  can  be  seen  more 
clearly  by  defining  a  parameter  as  follows: 


a 


(4.4) 


where  y^o  is  the  body  shape  given  by  the  first  order  outer  solution.  In 
the  range  of  values  of  n  considered  the  results  are: 


V  +  1 
2y  +  1 


5  n 


1 


a  =  0 


(4.5) 


Thus,  for  n  2  (y  +  l)/(2y  +  1),  the  entropy  layer  does  not  effect  the 
body  shape  obtained  from  the  first  order  outer  solution.  This  agrees  with 
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a  result  previously  obtained  by  Sychev  (Reference  3).  For  n  (y  +  1)/ 

(2y  +  1),  the  effect  of  the  entropy  layer  on  the  parameter  a  vanishes  for 
large  distances  downstream.  As  n->l/2,  the  effect  vanishes  more  slowly 
in  the  downstream  direction.  For  example,  withy  =  1.4: 


n 


n 


n 


=  3 

P 

II 

o 

4 

-0.520 

=  7 

Q  ~  X 

1 

-0.  143 

=  2 

a  ~  x 

(4.8) 


The  effect  of  nose  bluntness  on  the  inviscid  flow  field  can  be  illus¬ 
trated  by  considering  a  family  of  finite  length  bodies  of  the  same  thickness 
ratio. 


T 


<Vmax 

t 


(4.9) 


and  different  values  of  the  parameter  n.  A  family  of  such  bodies  with 
t  =  1/4  is  shown  in  Figure  2.  The  flow  field  profiles  at  the  base  station 
are  shown  in  Figures  3  and  4  plotted  versus  the  proportional  distance 
across  the  shock  layer: 


y* 


(4.10) 


The  existence  of  the  entropy  layer  is  most  clearly  illustrated  by  the 
velocity  and  entropy  function  profiles  of  Figure  4.  The  more  blunt  bodies 
or  lower  values  of  n  produce  thicker  entropy  layers  with  smaller  gradients 
at  the  body  surface. 
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21 


Figure  2.  Bodies  With  Power -Law  Shocks 


22 


23 


f).  CONCLUSIONS 


The  method  of  inner  and  outer  expansions  can  be  used  to  treat  the 
entropy  layer  near  the  body  for  blunt  nose  shpaes  with  power -law  shocks, 
y  ■ — xn,  for  1/2  5n  <1.  The  outer  solutions  near  the  shock  wave  are 
identical  to  the  blast  wave  solutions.  The  method  essentially  adds  a  con¬ 
stant  pressure  entropy  layer  to  the  blast  wave  solution  thus  providing  a 
solution  that  is  uniformly  valid  across  the  shock  layer. 

In  this  report  the  outer  solutions  were  expanded  about  the  body  to 
obtain  the  form  of  the  inner  expansions  for  1/2  Sn  *=  1 .  However,  as 
pointed  out  by  Lees  and  Kubota  (Reference  1),  for  n  «c  1/2,  the  outer 
solution  does  not  exist.  It  is  not  clear  how  the  method  of  inner  and  outer 
expansions  can  be  applied  to  the  case  n  <  1/2  or  if  such  an  analytic  solu¬ 
tion  can  be  obtained. 

The  inner  expansion  for  the  space  variable  y(x,  4<)  consistent  with 
the  first  order  outer  solutions  is  a  series  which  converges  very  slowly 
as  n— >1/2  from  above.  The  case  n  =  1/2  is  singular  and  requires  only 
one  term  in  the  expansion.  This  case  was  solved  previously  by  Yakura 
(Reference  4).  Also,  for  (y  +  l)/(2y  +  l)<n<l,  only  a  single  term  is 
required  in  the  expansion  and  the  inner  solution  is  easily  obtained.  In 
this  range  of  n,  the  blast  wave  profiles  are  modified  by  the  inner  solu¬ 
tion  but  the  body  shape  is  unchanged.  This  minimum  value  of  n  above 
which  the  blast  wave  body  is  correct  was  previously  presented  by  Sychev 
(Reference  3).  If  two  terms  are  retained  in  the  expansion  for  y(x,  ij>)  the 
inner  solution  can  be  extended  to  cover  the  range  (y  +  2)/(3y  +  2)  2  n 
c(y  +  l)/(2y  +  1).  However,  for  n  only  slightly  greater  than  1/2  the 
method  is  of  limited  usefulness  due  to  the  slow  convergence  of  the 
expansion. 
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